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In this paper we consider the issue of the Froissart bound on the high energy behaviour of total cross sections. 
This bound, originally derived using principles of analyticity of scattering amplitudes, is seen to be satisfied by 
all the available experimental data on total hadronic cross sections. At strong coupling, gauge/gravity duality 
has been used to provide some insights into this behaviour. In this work, we find the subleading terms to the 
so-derived Froissart bound from AdS/CFT. We find that a (In term is obtained, with a negative coefficient. 
We see that the fits to the currently available data confirm improvement in the fits due to the inclusion of such a 
term, with the appropriate sign. 


I. INTRODUCTION 

The high energy behaviour of the total cross sections for 
the scattering of any two high energy particles has been a 
subject of great theoretical interest over many decades, be¬ 
ginning from Heisenberg m. The various bounds one obtains 
on the rate of the rise of these cross sections are based just 
on the analyticity properties of scattering amplitudes and do 
not use any knowledge of the underlying dynamics of the in¬ 
teraction responsible for the scattering. The most important 
of all these has been the Froissart bound GllSl, which states 
that at high energies, the total cross section for a scattering 
process 1 + 2 —> / (/ = any final state) has an upper bound 
^ (In^ ^). Here y/s is the centre of mass energy and 
is an energy scale. Since the bound can be derived from very 
general physical arguments, such as the unitarity of the S ma¬ 
trix and certain analytical properties of the scattering ampli¬ 
tude, it has to be true in any quantum field theory. However, to 
derive this in the context of theories of strong interactions, like 
QCD, one requires to handle the theory in the non-perturbative 
regime. As a result, there exist only models and these usu¬ 
ally try to incorporate known properties of QCD, the mod¬ 
elling aspect involving assumptions and ansaetze about the 
non-perturbative regime iHlll. In fact, analyses of mini 
in the framework of the Bloch-Nordsiek improved eikonalised 
mini jet model Ha indicate a direct relationship between the 
Froissart bound at high energy and the dynamics of ultra soft 
gluons, i.e., the behaviour of QCD in the far infrared. An anal¬ 
ysis of the high energy behaviour of the available data on to¬ 
tal hadronic cross sections, while being completely consistent 
with the bound, seems to need in the fits the presence of a 
subleading (In j-) term to the Froissart bound. Various QCD 
(inspired) modefs use such a subleading term, though there 
is no theory with a proper explanation of this subleading be¬ 
haviour miiiHibi. 

In this article, we attempt to seek a theoretical understand¬ 
ing of this subleading behaviour shown by the data on total 
cross sections at high energies, using the AdS/CFT correspon¬ 
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dence ini. It is known that the leading term (Froissart bound) 
can be derived using AdS/CFT ifTSltSOll . The model based 
on lEl describes high energy scattering in a large-W,Iarge- 
A gauge theory with broken conformal symmetry which is 
achieved by putting a cut-off in the infra-red (IR). Here N is 
the number of colours and A is the’t Hooft coupling. When a 
point mass m is placed on the IR boundary, the perturbations 
in the AdS space can be big enough to create a black hole. The 
geometrical cross section of such a black hole, whose radius 
is equal to the AdS radius, is the gravity dual of the maximum 
possible scattering cross section in the field theory. In order 
to have a more realistic model, which takes into account finite 
coupling corrections (in 1/N and 1/A), we can incorporate 
higher curvature corrections in the dual gravity description. If 
we are working perturbatively in the couplings, then the lead¬ 
ing correction of this type will be at four derivative order ll22l . 
Such terms can be redefined into a single Gauss-Bonnet term 
whose coefficient describes a 1 /N correction - see eg. Il2^ for 
the dictionary relating the field theory and gravity in such a 
theory. As we will argue below, this result will in fact hold for 
any higher curvature gravity dual based on the results in ll24l . 
In ifT^I^ it is argued that the corrections due to string modes 
are exponentially small, so here we will ignore them from the 
onset. However it is important to take into account 1 /N and 
1/A corrections to ensure that the behaviour of the subleading 
terms arising from Einstein gravity as additional contributions 
to the Froissart bound do not get significantly modified. 

In this paper, we will consider subleading corrections to the 
Froissart bound that arise from such an analysis using the re¬ 
sults of M- The incorporation of the higher curvature terms 
in the dual gravity will enable us to see what the finite cou¬ 
pling effects are. We show that the subleading term in that 
cross section is indeed (In ^), as seems to be required by the 
fits. Our AdS/CFT analysis predicts that the coefficient of this 
subleading term should be negative. Further, if we consider 
the subleading term coming in the dual impact parameter from 
Einstein gravity, in the dual field theory there is an additional 
correction to the Froissart bound of the form (In —) (In In —). 
A recent work II25I also found the same subleading term as a 
correction to the maximal impact parameter. At high energies, 
we expect that the general structure of the subleading terms 
that we find from AdS/CFT will be useful in getting better fits 
to experimental data. This is precisely what we find. 

Finally, we use all available data (including FHC and cos- 
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mic ray data) to conclude that the (In term significantly 
improves the fit (the F test is used). Also, the data always 
chooses the coefficient of the (In ^) term to be negative, con¬ 
firming our prediction from holography. 


The stress-energy tensor Tab is generated by any matter 
present and is given by 


Tab 



5Cm 

ggab 



(7) 


II. ADS/CFT DERIVATION OF THE EROISSART BOUND 


In this section, we will review and redo the analysis in ifTSl 
relating the gravity calculations to the Froissart bound in order 
to find the effects of finite coupling. A more careful analysis 
of this can be performed using the results in ll20l . Since the 
functional form of the transcendental equation which gives the 
maximal impact parameter is similar to what we will use, we 
do not expect that the conclusions we arrive at would change 
significantly in this case. 


B. Solution as background and perturbations 

To solve the equations of motion (|^, we introduce the met¬ 
ric ansatz 

3ab — 3ab “t“ Flab^ (8) 

where Hab are the perturbations on an AdSs background met¬ 
ric gab given by 


A. Action and equations of motion 



= diag{-l,1,1,1). (9) 


The most general quadratic gravity action in 5 dimensions 
is given by 

S = Seh + Sgb, (1) 

where Seh is the standard Einstein-Hilbert action with a cos¬ 
mological constant A 


Seh = 




( 2 ) 


and Sgb is the Gauss-Bonnet higher derivative correction to 
it given by 


z is the radial coordinate of the AdSs and z = L is its IR 
boundary: 0 < z < L. We use the convention that indices 
with greek alphabets (g, v ,range from 0 to 3 and are raised 
(lowered) with the metric Indices with latin alpha¬ 

bets (a, b,...) include the radial coordinate z as well and are 
raised (lowered) with 5 “^ {gab)- The background metric gab 
satisfies the Einstein’s field equations with a negative cosmo¬ 
logical constant A = 

1 6 

Rab l^Ttgab 'J^do.b — 0 - ( 10 ) 

We redefine the perturbations to a convenient form 


^gb = IJ d^x^g - m-'^Rab + le^^’^Rabcd) , (3) 

with K = 7 being a perturbative coupling. The 

5-dimensional gravitational constant GT’'> is related to the 
Planck mass Mp as 



and assume that the perturbations are small enough to set 

{hab)'" « 0 . ( 12 ) 


G(5) 


1 

= Mf- 


(4) 


gab is the metric and p = det(pa 6 ). i? is the Ricci scalar, is 
the 5-dimensional Ricci tensor and Rabat is the 5-dimensional 
Riemann tensor. The term Cm is the contribution due to any 
matter fields in the theory. 

Extremizing this action gives the 5-dimensional equations 
of motion 


We further assume that 


hzz = 0 = hp,:,. (13) 

Therefore we are left with ten nonzero components of hab- 
But there are 15 equations of motion and hence we have the 
freedom to choose the traceless-transverse gauge 

= 0, h = = 0, (14) 


G ab F gab A -f ^Rab — blTab, (5) 

where Gab = Rab — \R9ab is the Einstein tensor and Bab is 
defined as 


where is the covariant derivative in the AdSs background. 

The equations of motion for the perturbation can be derived 
from 0 using ( [T 0 | . In the traceless-transverse gauge, the lin¬ 
earised equations of motion for hpi, are 


Bab 


-^gab[RcdefR‘^^^^-^RcdR’^‘^ + R")+2RRab = ^ = ( 15 ) 

+ 2RacdeRb^‘^^ - 'iR^abRd" " ^RaA"' 

( 6 ) where □ = g‘^’’'Va^b, using the unperturbed metric. 
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Thus the effect of the 1/7V correction will be to simply 
rescale the Newton constant. AdS/CFT arguments suggest that 
7 > 0 and so, C < 1 although 7 < 0 cannot be ruled out. 
Moreover, in this theory, the first 1/A correction arises at eight 
derivative order through R'^ corrections (see eg. l 23 ), the re¬ 
sult of which from the gravity side will be further subleading. 
Since we are treating the curvature corrections perturbatively, 
we can in fact give a stronger argument for the generality of 
our result based on Il24l . Notice that in the arguments above, 
we need to find the linearised equations around AdS. In ll24ll it 
was shown that a very general higher curvature action can be 
rewritten in the form of a background field expansion around 
AdS. When this is done, the linearised equations of motion 
will only be sensitive to the quadratic curvature terms in the 
background expanded action. Since we are working pertur- 
batively, we can use field redefinitions to get back the Gauss- 
Bonnet gravity. Hence our claim above will work quite gen¬ 
erally with the coefficient C related to the two point function 
of the stress tensors in the field theory. In the general case, C 
can be either > 1 or < 1 so long as it is positive, which is 
demanded by unitarity of the field theory. 

The rest of the analysis thus parallels M, which we will 
review below in order to understand the systematics of the sub¬ 
leading terms in the Froissart bound. 


III. SOLUTION TO THE PERTURBATIONS 


We are interested in solving for the linearised perturbations 


which satisfy (151 with a point mass m on the IR bound¬ 


ary z = L. This point mass is a source and its contribution 
in the action is accounted for by adding a suitable Cm- The 
addition of this mass leads to a stress tensor with only one 
non-zero component 


Tab = 


non-zero 


if a = & = 0, 
otherwise. 


( 16 ) 


But such a Tab is not traceless 

g^^^Tab ^0 (17) 

and hence is not compatible with the equation of motion ra¬ 
in order to make it traceless, we also add an incompressible 
gas on the brane, which only generates a pressure Tn and no 
shear stress (T^- = 0 for i ^ j). The resultant traceless stress- 
energy tensor is 

Tab= mS^(x)S{z-L){6°S^ + 6lSl) . (18) 

With this source, the only non-trivial equations of motion are 

Dhii =S^(x)S{z-L), j = 0,1. (19) 

But before we solve 0, we need to specify the boundary 
conditions. As the /i^^’s are gravitational fields, they must 
satisfy Neumann boundary conditions in the IR brane 

( 20 ) 


where is the unit vector outward, normal to the IR bound¬ 
ary. 

The boundary problem ( [T^ can be solved using the scalar 
Green function (X, X') satisfying 


(A, X') = (x - x') S {z - z '), (21) 


with coordinates X = (z,x). The boundary condition (201 
translates to 


a.Alf ix,x') 


= 0. 


2=L 


( 22 ) 


f 5) 

The perturbations ha can be obtained from ^ as follows: 
hu = ~ I S^{x\z'-L), (23) 


where 

S'^{x', z’ -L) = - L). (24) 

Even though is a second order partial differential equa¬ 
tion, in the Fourier space it reduces to a second order ODE in 
the variable z. The Eourier transform in Minkowski 4 dimen¬ 
sions defined as 

A(f (25) 

J y 2i7r J 

must satisfy 

{z'^dl-3zd^ Tq'^z'^) Ap = 5{z-z’), (26) 

with = —p^. Under the redefinition 

equation ( |26l l becomes 

{z'^dl + zdz + q^z^ — 4 ) Ap = L—S (z — z'). (28) 

z' 

Eor z ^ z', the above equation admits as its two independent 
solutions the Bessel functions J 2 (gz) and Y 2 {qz): 


Ap (z,z') 


A+ (z, z') if z > z' 
A_ (z, z') if z < z', 


(29) 


where 


A± (z, z') = A± (z') J 2 (qz) + B± (z') Y 2 (qz). (30) 


A± and B± are functions of z' and can be determined from 
the boundary and matching conditions. 

In the region z > z', the boundary condition (22 1 sets 


A+ = -B+ 


YijqL) 

MqL)- 


( 31 ) 
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In the region z < z', demanding that be regular at z = 0 
yields 


= 0 . 


(32) 


Equation (281 implies Ap is continuous across z = z' 


A+(z,z') =A_(z,z') 


(33) 


but its derivative in z is not. The value of the discontinuity in 
9zAp is easily found by integrating (28 1 across z = z': 


dz 


A+ {z,z') - A_ (z,z') 


= (34) 


The matching conditions (33 i and ([34|i yield 


A _ ttL Ji(qL)Y2(qz')-J'i(qz')Y-i(qL) 
- “ 2 Ji(9L) 

B+ = ^J2{qz'). 


(35) 


Consequently, we can express Ap as 


^ ttL / z<z> n 2 J2(gz<) 

2 \ J Ji{qL) 

[Ji{qL)Y2{qz>) - Yi{qL)J2iqz>)] , (36) 


where z< (z>) denotes the lesser (greater) of z and z'. 

The Green function A(^^ {X,X') can be obtained by the 
inverse Fourier transform 


Aif iX,X') 


ttL / z<z> n2 r d‘^p J2iqz<) 

2 \ ) J (2^)4® J^{qL) 

[MqL)Y2{qz^) - Yi(gL)J2(9^>)] ■ (37) 


In the case where one of the arguments of A ^^' is on the IR 
brane, at z' = L, one can use the Wronskian of the Bessel 
functions 


r„(gL)J„+i(gL) - J„(gL)y„+i(gL) = (38) 

valid for all n G N to reduce ([3^ to 


Using (23 1 and integrating over x' and we get 

■Z\^ f (Pp „ip.g J 2 MZ) 






(40) 


Although the integral in ( |40] l is difficult to evaluate explic¬ 
itly, we are only interested in the long-distance z « L limit 
where it simplifies. In this case, the integral is dominated by 
the region of small \(f\z and we can replace the Bessel function 
J 2 {\<f\z) by a small argument expansion. We find 


hu (A) 


mn ( 


d^P 

(2^)3 Ji(lqlL)- 


(41) 


In the spherical coordinates 

d^p = 27r|p|^c?|p|d (cos 9 ), 
integration over the angular dependence yields 


( 42 ) 


(A) 

rriK 
327r2C' V L 


1 


poo pi\p\r _ p-i\p\r 


where r = |a:|. 

The denominator of (|4^ has zeros where \q]L satisfies 


Ji ML) = 0. 


(44) 


There are infinite number of such zeros and hence the inte¬ 
grand has infinite number of first order poles. These poles are 
located on the positive real axis of the -plane: 

k1poie = ^, ji.feeK+, fc=l,2,... (45) 

It is easy to see that in the long-distance (large r) limit, the 


contributions from the higher poles in the integral (43 1 are ex 


ponentially suppressed and hence can be ignored. Thus, it suf¬ 
fices to consider the first pole at \q\L = ji,i, but we will keep 
the second pole at \q]L = ji 2 as well to explicitly check the 
above claim. 

The denominator of ( |43] l can be expanded in the small 
neighbourhood of the poles, where the integrand is expected 
to contribute significantly. Using the asymptotic behaviour of 
the Bessel functions for a small argument, we can write 


hm Ji(|^L)« ---. 

z 


(46) 


This property can be used to evaluate the residues at the poles 
jip and ji_ 2 . The integral can thus be evaluated to 


16nCL 


(I) 


4 p 


i-f 


Jl,2 

Jl.l 


r 


(47) 


The second term in the parenthesis is the contribution from 
the second pole ji 2 - As previously suggested, this term is 
exponentially suppressed as compared to the leading term in 
jip and hence can be safely ignored in the large r limit. 


IV. GAUGE/GRAVITY DUALITY: UPPER BOUND TO THE 
SCATTERING CROSS SECTION 

A sufficiently large mass m in ( |47) l will make the pertur¬ 
bation big and when hu ~ 1 we can expect a black hole to 
be formed m. The causal horizon of the black hole lies at 
r = th, where iJoo = 1 (where goo vanishes). 
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So, in the IR boundary at z = L, an estimate for rn can be 
obtained from (|47]i as 


16ttCL 




rH 


oc p, 


(48) 


with p a constant of order 1: p ~ 0(1). Identifying the black 
hole energy with its mass, E (y: m, we can re-express ([48]l as 


rn 


Ml 


E 


th 


— In — -In^ -InC 


En 


(49) 


In the high energy (large s) limit (53 1 reduces to 
'1 


< 


Mf 


lir -h /3 In-(In In —) In — 

4 So So So So 


(54) 


This is our main result, which we will use in the fits to exper¬ 
imental data below. The sign of the coefficient /3 needs to be 
determined next. Here, 



(55) 


where Eq = has the dimension of energy. Mi = 

is the mass of the lightest Kaluza-Klein mode of the graviton 
and (7 > 0(1) is a constant related to C which absorbs the 


ambiguity in the RHS of (48 1 . 


To solve the transcendental equation (49 1 , we insert the 
ansatz 


I ( E ( E 


(50) 


with Cl, C 2 constants to be determined. In the long-distance, 
high-energy (r, E large) limit, ( |50| l satisfies the equation with 

1 1 
^^~~Wi “ MiL' 

The geometric cross section of the black hole is then given by 
(^BH = ujirrjj, (51) 


where oj is some constant different from unity which takes 
into account the higher curvature effects through the Wald for¬ 
mula. We can absorb this factor into 1 /Mf, which sits outside 
our formulae - this does not play any significant role in our 
analysis of the general structure. 

This gives an estimation of the maximal scattering cross 
section in the gauge theory ca. The cross section in the gauge 
theory is bounded from above by the geometric cross section 
of the black hole: 


O’ < CTbh- 


(52) 


Since we want to use the standard notation s = E'^ and nor¬ 
malise using So = 2 mp, with uip being the proton mass, we 
can set {E/Eq)'^ = (s/so)e^ to get 


< —K ( - In-h- — InC — ln( 

~ ]\/f2 I o ^ ' o V 


Ml 


So 


^ -[ln-+C]) 


2MiL sq 


(53) 


The first zero of (44 1 is at ji i « 3.8, and so MiL « 3.8. 


V. CORRECTIONS TO THE EROISSART BOUND 


In gauge/gravity duality, M^ = 2iV^/ (ttL^) where N is large 
and thus mp/Mp <C 1. Further MiL « 3.8 with C > 0(1), 
as mentioned earlier. This means that the argument of the log¬ 
arithm is ~ rripL'^/N'^C'^. Hence we expect /3 < 0, since 
the argument of the logarithm is <C 1 for large N. Thus the 
main conclusions from our analysis are that the structure of 
the subleading terms which arise from the Einstein dual are 
not altered at finite 1/A^, 1/A and that the subleading term 
(In j^) should come with a negative coefficient. Notice that 
in order to reach this conclusion we just needed sq to be less 
than the scale Eq arising naturally in the AdS/CFT calculation 
and the finite coupling effects (in 1/iV and 1/A) to contribute 
perturbatively (so that C does not become drastically small 
for instance). We have two predictions from AdS/CFT. First 
that the sign of the (In ^) is negative and second that the ra¬ 
tio 1281 between the coefficients of the (In^ ^) term to the 
(In ^)(lnln ^) term is —1/4. We will now confront experi¬ 
mental data with these predictions. 


VI. EITS TO EXPERIMENTAL DATA 


The experimentally measured total hadronic cross section 
for pp and pp scattering events shows an initial de¬ 

crease and a later rise with the centre of mass energy An 
understanding of this form is beyond the scope of perturbative 
QCD. The theory of analyticity of the S matrix demands that 
the hadronic amplitudes be analytic functions in the complex 
angular momentum J. In this plane, the singularities do not 
depend on the scattering hadrons. Unitarity of partial waves 
and boundedness of the absorptive part within the Lehman 
ellipse lead to the Froissart (upper) bound on Also, 

Regge trajectory exchanges and a Pomeron exchange seem to 
play an important role in the description of 

The behaviour of the total cross section was first success¬ 
fully parametrised by Donnachie and Landschoff (DL) ll29l 
[30l. Using Regge-Pomeron theory, they proposed a sum of 
two powers to describe the then available data: 


As we pointed out in the introduction, the (In ^) term has 
been widely used in fits to experimental data. Recent develop¬ 
ments in the theory suggest that a (In ^)(lnln ^)-like term 
might be necessary in addition at higher energies yTl- 


= + (56) 

where X, T, e, and 77 are parameters and X should be equal for 
cr//) and But as new data in a higher energy regime be¬ 
came available, the DL parameterisation became insufficient. 
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^2 

^3 


a (mb) 

36.980 

30.692 

37.398 

31.786 

b (mb) 

38.3 ± 1.6 

50.0 ± 1.5 

36.87 ±0.70 

45.38 ±0.20 

c (mb) 

-1.42 ±0.12 

0 

-1.44 ±0.11 

0 

d (mb) 

0.043 ± 0.042 

0.153 ±0.048 

0 

0 

e (mb) 

0.276 ±0.011 

0.1687 ±0.0075 

0.2844 ±0.0076 

0.1922 ±0.0017 

/(mb) 

-21.5 ± 1.0 

-18.7 ± 1.0 

-20.95 ±0.85 

-16.78 ±0.71 

a 

0.490 ±0.012 

0.525 ±0.014 

0.496 ±0.010 

0.550 ±0.011 


270.535 

276.018 

270.606 

391.754 

Ndf 

342 

343 

343 

344 

Ndf 

0.791 

0.805 

0.789 

1.139 


TABLE I: Values and errors of the fit parameters in 1 57 i-i 60 1 . The table includes the corresponding values of x 




Arguments from crossing symmetry, analyticity and unitar- 
ity have since been used to modify the DL parameterisation 
and provide a satisfactory description of the data (for exam¬ 
ple, see EiEiini). An understanding of pp and pp scattering 
is not only intrinsically interesting, but also leads to impor¬ 
tant insights into 7p and 77 (7 = photon) scattering processes 

(eg. IHEII). 

It is obvious that (j5^ will be a better fit to experimental 
data than just the (In^ term, because of the bigger parame¬ 

ter space available. But the question one should ask is whether 
the improvement to the fit is just because there are more pa¬ 
rameters in the theory or whether the inclusion of the correc¬ 
tion terms is really necessary to improve the fit. The answer 
to this question can be given by figuring out whether the terms 
are statistically significant in the F test. 


A. Data fitting 

A good parametrisation of the data separates the cross sec¬ 
tion into two parts: a piece which saturates to a constant value 
and a piece which accounts for the rise ifTSll . The first piece 
contains a constant term due to the Pomeron and a Regge term 
decreasing as 1/ y/s. Two such Regge terms are required, rep¬ 
resenting the C-even and C-odd exchanges. It is well-known 
that oc 13^ . The second piece comes from 

a triple pole at J = 1, which produces (In^ ^) and (In ^ ) 
terms in the total cross section m. At large energies, the 
Regge terms do not contribute. The above mentioned logarith¬ 
mic piece certainly dominates asymptotically, but plays an im¬ 
portant role at all energy values. This is due to the fact that the 
fitted functional form has to describe the early fall of the cross 
section with the energy as well as the common (for pp and pp) 
constant value from where the high energy rise begins. Re¬ 
cently, the determination of the energy scale of the logarithmic 
terms has been shown to further contribute a (In ^) (In In ^) 
term to the total hadronic cross section fT7\ '. Consequently, 


^ This term is also relevant in the case of inelastic cross sections (H 


we have simultaneously fitted experimental data points 

to the following functional form: 

af = a -\—^ -I- clny -I- din j/(ln(lny)) + e\v? y ± 

(57) 

where the upper (lower) sign refers to pp (pp), y = 

and {a, b, c, d, e, /, a] are fitting parameters A single set 
of these parameters accounts for both pp and pp data. 


It is important to note that the term ±fy°‘~^ in (57 1 plays a 
fundamental role in the simultaneous description of pp and pp 
scattering. This term accounts for the difference between pp 
and pp scattering in the low energy region while ensuring that 
a unique set of parameters is used in the high energy region, 
where there is little difference between the two sets of data ^. 

We also consider simpler fits by constraining the value of 
some of these fitting parameters: 


^2 


(7 q - (7 


s.t. d = 0 , 


at' = at s.t. c, d = 0 . 


(58) 

(59) 

(60) 


The analyticity of these u^’s enforces the constraint m 
a = 48.58 - 0.35166 - 2.091c -f 2.715d - 4.371e. (61) 


The fit ( [57| ) takes the functional form used by m and adds 
the term dln 2 /(ln(ln?/)) to it. The fits (58 1 - ( 6 O 1 are required 


^ The parameter a is sometimes set to a = 1/2 (for example in (H.). Here 
we do not assume a = 1/2, but rather we show from the fits that a ^ 1/2 
(see TABLEig. 

^ Independentnts to pp and pp data without lead to different val¬ 

ues of the parameters c, d, e relevant to the description of the data at high 
energies in the two cases. Since we want to describe both pp and pp data 
with the same high energy parameters, we do not perform such separate 
fits. It should be stated however, that the conclusion about the sign of the 
(In ^) term that we draw in the end is unchanged even then. 
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LogjuVs" (GeV) 

FIG. 1: (Colour online.) Fit results to experimental values of and The magenta solid, orange dot-dashed, blue dashed and black 
dotted curves are the - j60| l fits to the pp (green circles) and pp (red squares) data points, respectively. The data are from CDF, E710, E811, 
UAl, UA4, UA5 experiments f35P^ . The pp data points also include results from the LHC (at y/s = 7, 8 TeV) l43ti45l and cosmic-ray 
data 1461 . 


to study whether correction terms found from gauge/gravity 
duality in the previous section (see equation ( |54) i) lead to a 
better description of the data. 

Setting v/Snjin = 6 GeV allows us to exploit the rich sample 
of low energy data just above the resonance region. 

The goodness of the fits has been estimated with the chi- 
squared value per degree of freedom: 


_ 


X 


Ndf Nd — Np 


(62) 


where Nd is the number of data points considered and Np is 
the number of parameters in the fit. The chi-squared value 
has been taken to be 


/ (fit) (exp)' 


(63) 


where is the total cross section given by the fit, is 

the experimentally measured total cross section and 5ai is the 
error in the experimental total cross section. If the chi-squared 
value per degree of freedom is of order 1 


X 


2 


Ndf 


6 ^( 1 ), 


(64) 


then we consider the fit to be good. 

The results of the fits ([57|-([60| constrained by ( [6T] l are tab¬ 
ulated in TABLE 1^. All fits are good, since they satisfy ( [64| l. 
Note the small errors in the fit values of c. Thus c is indeed 
negative to a very high degree of significance. 

The pp and pp cross sections derived from the parameters 
of TABLElJare shown in EIG.[T]as a function of the centre of 


mass energy y(s. The pp data points (red squares) include 
crffj results from the experiments CDF, E710, E811, UAl, 
UA4 and UA5 ll35H42ll . The pp data points (green circles) in¬ 
clude crfot results from LHC (at = 7,8 TeV) B3lj45l and 
cosmic-ray data ll46l . The magenta solid, orange dot-dashed, 
blue dashed and black dotted curves are the , cr^, crj 

fits to the pp and pp data, respectively. 


B. Comparison of the fits 


In this subsection, we study whether and which one of the 
fits (57 1 -([ 60 | provides a statistically significantly better de¬ 
scription of the experimental data. 

The fit ( [60| l is the simplest of all: it can be obtained from 
by setting c = 0 or d = 0 or both. Hence, it can be 
thought of as a special case of any of Similarly, ( |58| ) 

and (59 1 can be obtained from (|5^ by setting c = 0, d = 0, 
respectively. 


For this reason, we can consider pairs of such fitting mod¬ 
els as nested models. One can compare a pair of models which 
are nested by taking one of them to be the null hypothesis (the 
model with fewer parameters) and the other one as the alter¬ 
native hypothesis (the model with more parameters). It is ob¬ 
vious that the alternative hypothesis will almost always have a 
2 

smaller —, due to the bigger parameter space available. In 

2 

order to determine if such an improvement in — is signif- 
leant, the F test needs to be used (for a review of the F test, 
see l47ll ). 

In a U test, the and Np,p values of the null and altema- 
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tive hypotheses are used to compute the F ratio as 


F = 


^^2(null) _ ^2(alt)^ 


(65) 


If F « 1, then the improvement in provided by the 
alternative hypothesis should not be regarded as significant. 
Rather, the null hypothesis should be considered sufficient to 
account for the data. The converse is true if F > 1. 

Since the F distribution is known, using the F ratio and the 
quantities 


Nn = ^ . Nd = N^df > ( 66 ) 

a P value can be obtained. Such P value quantifies the prob¬ 
ability with which one expects a particular F ratio value (or 
higher), provided the null hypothesis is correct. If F < /3, 
with P the significance level chosen, then the alternative hy¬ 
pothesis should be considered to fit the data significantly better 
than the null hypothesis, and vice-versa. 

TABLE |I^ contains the relevant quantities in the computa¬ 
tion of P for the comparison of all pairs of nested models. It 
is easy to see that must be rejected in favour of any other 
model, with a significance level /? = 0.0001 when compared 
to af and and with a significance level /3 = 0.002 when 
compared to ■ <jf is better than tr^ with a significance level 
P = 0.0001. There is no statistically significant difference be¬ 
tween crj*^ and ■ 


Null 

Alt 

F 

Nm 

Nd 

P 

^2 


140.7780 

1 

342 

< 0.0001 



1.0522 

1 

342 

0.3057 

^4 


77.3169 

2 

342 

< 0.0001 

0-4 


9.8444 

1 

343 

0.0019 

^4 

^3 

153.5580 

1 

343 

< 0.0001 


TABLE II: F and P values for the comparison of all possible pairs 
of nested models. The null hypothesis is the simpler model. The 
alternative hypothesis is the model containing more parameters. 

Hence, using the logic of the F test mentioned above, we 
conclude that ( |59| l provides the best description to the exper¬ 
imental data among the functional forms considered here. A 
functional form including a logarithmic term (In ^) produces 
a statistically significantly better description of the data than 
just a Froissart-like bound term alone. The coefficient of this 
subleading term is negative. The inclusion of a further sub¬ 
leading term of the form (ln^)(lnln^) does not signifi¬ 
cantly improve the fit. Hence, as we previously pointed out, 
the inclusion of more parameters does not necessarily improve 
the data description significantly, as it is here the case. 

It can be argued that the (In ^)(lnln correction term 
may become distinguishable at higher energies, where it con¬ 
tributes most. Indeed, in the energy range where data are avail¬ 
able, this term behaves almost as a constant and is sublead¬ 
ing to the (In term for the values of c, d found in (see 


TABLE|^. More and better data from air shower experiments 
may be the only possibility of ever addressing the issue from 
the point of view of achieving any further discrimination. 


VII. DISCUSSIONS 

The high energy behaviour of the total hadronic cross sec¬ 
tions was first analysed in m, where the bound on the rise 
of the cross section was shown to arise from the finite range 
of the interaction. In the case of QCD, this short range arises 
from the properties of the theory in the infrared. This indi¬ 
cates that it is the behaviour of the theory of strong interac¬ 
tions in the IR that decides the high energy behaviour of the 
cross section. Arguments from analyticity and unitarity lead 
to the Froissart (upper) bound on such total cross sections. In a 
Bloch-Nordsiek improved eikonalised mini jet model lfT3l . the 
Froissart bound seems to be directly related to the behaviour 
of QCD in the far IR. Phenomenological fits usually include a 
subleading logarithmic term as well. A logarithmic subleading 
term is often used to fit data, but lacks a theoretical explana¬ 
tion. It is just a phenomenological observation. 

Ideally, one would like to understand such a subleading 
term in the context of QCD. The understanding of the be¬ 
haviour of total hadronic cross sections lies in the realm of 
non-perturbative QCD. Although there has been a lot of in¬ 
terest in the subject and many efforts have been made, the 
tools required for handling non-perturbative QCD are still be¬ 
ing developed. However, using the AdS/CFT correspondence, 
one can map a problem in a strongly coupled gauge theory 
to a problem in a weakly coupled gravity theory. Hence, we 
have mapped the high energy behaviour of the total hadronic 
cross section to a gravity toy model and investigated what the 
correction terms to the Froissart bound are in this simplified 
scenario. 

We have shown that the extension of the holographic argu¬ 
ments in ifTSl [191 generates a subleading (In term to the 
Froissart bound in the dual field theory. The duality predicts 
that the coefficient of such a subleading term is negative if 
the finite coupling corrections are small. In other words, the 
prediction from an Einstein dual together with a perturbative 
1/A^, 1/A correction through higher curvature corrections is 
that the coefficient of (In is negative. Therefore (since it 
reduces the upper bound), it is an improvement to the Froissart 
bound. 

Further, we have used all available data to demonstrate 
that such a (In term indeed significantly improves the fits. 
Also, the prediction about the negative coefficient for the sub¬ 
leading term is confirmed from the fits. 

Presumably, the (In^) term reflects the nature of non- 
perturbative QCD. While one of course cannot claim that the 
holographic theory is modelling QCD, it is quite remarkable 
that the general structure arising from holography for the sub¬ 
leading terms in the Froissart bound as well as the sign of the 
first subleading term agree so well with data. Hence, it is pos¬ 
sible that the gravity toy model here considered captures cer¬ 
tain important aspects of non-perturbative QCD. 

While the argument from AdS/CFT leading to the negative 
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coefficient for the (In —) term seems reasonable, another pre- 
diction from our AdS/CFT analysis is that the relative coeffi¬ 
cient between the (In^ —) term and the (In — )(lnln —) term 
is -1/4. It would be interesting if this was home out by the 
data, but it seems unlikely to happen in the near future. 
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